Coupling remote sensing and eDNA to monitor environmental impact: A pilot to quantify the environmental benefits of sustainable agriculture in the Brazilian Amazon.

Monitoring is essential to ensure that environmental goals are being achieved, including those of sustainable agriculture. Growing interest in environmental monitoring provides an opportunity to improve monitoring practices. Approaches that directly monitor land cover change and biodiversity annually by coupling the wall-to-wall coverage from remote sensing and the site-specific community composition from environmental DNA (eDNA) can provide timely, relevant results for parties interested in the success of sustainable agricultural practices. To ensure that the measured impacts are due to the environmental projects and not exogenous factors, sites where projects have been implemented should be benchmarked against counterfactuals (no project) and control (natural habitat) sites. Results can then be used to calculate diverse sets of indicators customized to monitor different projects. Here, we report on our experience developing and applying one such approach to assess the impact of shaded cocoa projects implemented by the Instituto de Manejo e Certificação Florestal e Agrícola (IMAFLORA) near São Félix do Xingu, in Pará, Brazil. We used the Continuous Degradation Detection (CODED) and LandTrendr algorithms to create a remote sensing-based assessment of forest disturbance and regeneration, estimate carbon sequestration, and changes in essential habitats. We coupled these remote sensing methods with eDNA analyses using arthropod-targeted primers by collecting soil samples from intervention and counterfactual pasture field sites and a control secondary forest. We used a custom set of indicators from the pilot application of a coupled monitoring framework called TerraBio. Our results suggest that, due to IMAFLORA's shaded cocoa projects, over 400 acres were restored in the intervention area and the community composition of arthropods in shaded cocoa is closer to second-growth forests than that of pastures. In reviewing the coupled approach, we found multiple aspects worked well, and we conclude by presenting multiple lessons learned.

Seasonal development of a tidal mixing front drives shifts in community structure and diversity of bacterioplankton.

Bacterioplankton underpin biogeochemical cycles and an improved understanding of the patterns and drivers of variability in their distribution is needed to determine their wider functioning and importance. Sharp environmental gradients and dispersal barriers associated with ocean fronts are emerging as key determinants of bacterioplankton biodiversity patterns. We examined how the development of the Celtic Sea Front (CF), a tidal mixing front on the Northwest European Shelf affects bacterioplankton communities. We performed 16S-rRNA metabarcoding on 60 seawater samples collected from three depths (surface, 20 m and seafloor), across two research cruises (May and September 2018), encompassing the intra-annual range of the CF intensity. Communities above the thermocline of stratified frontal waters were clearly differentiated and less diverse than those below the thermocline and communities in the well-mixed waters of the Irish Sea. This effect was much more pronounced in September, when the CF was at its peak intensity. The stratified zone likely represents a stressful environment for bacterioplankton due to a combination of high temperatures and low nutrients, which fewer taxa can tolerate. Much of the observed variation was driven by Synechococcus spp. (cyanobacteria), which were more abundant within the stratified zone and are known to thrive in warm oligotrophic waters. Synechococcus spp. are key contributors to global primary productivity and carbon cycling and, as such, variability driven by the CF is likely to influence regional biogeochemical processes. However, further studies are required to explicitly link shifts in community structure to function and quantify their wider importance to pelagic ecosystems.

Universal closed-tube barcoding for monitoring the shark and ray trade in megadiverse conservation hotspots.

Trade restrictions for endangered elasmobranch species exist to disincentivise their exploitation and curb their declines. However, trade monitoring is challenging due to product variety and the complexity of import/export routes. We investigate the use of a portable, universal, DNA-based tool which would greatly facilitate in-situ monitoring. We collected shark and ray samples across the Island of Java, Indonesia, and selected 28 commonly encountered species (including 22 CITES-listed species) to test a recently developed real-time PCR single-assay originally developed for screening bony fish. In the absence of a bespoke elasmobranch identification online platform in the original FASTFISH-ID model, we employed a deep learning algorithm to recognize species based on DNA melt-curve signatures. By combining visual and machine-learning assignment methods, we distinguished 25/28 species, 20 of which were CITES-listed. With further refinement, this method can improve monitoring of the elasmobranch trade worldwide, without a lab or species-specific assays.

Resource competition drives an invasion-replacement event among shrew species on an island.

Invasive mammals are responsible for the majority of native species extinctions on islands. While most of these extinction events will be due to novel interactions between species (e.g. exotic predators and naive prey), it is more unusual to find incidences where a newly invasive species causes the decline/extinction of a native species on an island when they normally coexist elsewhere in their overlapping mainland ranges. We investigated if resource competition between two insectivorous small mammals was playing a significant role in the rapid replacement of the native pygmy shrew Sorex minutus in the presence of the recently invading greater white-toothed shrew Crocidura russula on the island of Ireland. We used DNA metabarcoding of gut contents from >300 individuals of both species to determine each species' diet and measured the body size (weight and length) during different stages of the invasion in Ireland (before, during and after the species come into contact with one another) and on a French island where both species have long coexisted (acting as a natural 'control' site). Dietary composition, niche width and overlap and body size were compared in these different stages. The body size of the invasive C. russula and composition of its diet changes between when it first invades an area and after it becomes established. During the initial stages of the invasion, individual shrews are larger and consume larger sized invertebrate prey species. During later stages of the invasion, C. russula switches to consuming smaller prey taxa that are more essential for the native species. As a result, the level of interspecific dietary overlap increases from between 11% and 14% when they first come into contact with each other to between 39% and 46% after the invasion. Here we show that an invasive species can quickly alter its dietary niche in a new environment, ultimately causing the replacement of a native species. In addition, the invasive shrew could also be potentially exhausting local resources of larger invertebrate species. These subsequent changes in terrestrial invertebrate communities could have severe impacts further downstream on ecosystem functioning and services.

Molecular phylogeny and systematics of bald uakaris, genus Cacajao (Primates: Pitheciidae), with the description of a new species.

Bald uakaris, genus Cacajao, are Amazonian primates currently classified as one species and four subspecies based on the patterns of pelage coloration. In this study, we test if their current taxonomy is represented by the phylogenetic relationship of the main lineages retrieved from molecular data. We included, for the first time, all bald uakari taxa in a mitochondrial (cytochrome b) and genome-wide (ddRAD) phylogenetic analyses. We also examined the pattern of pelage colouration in specimens from zoological collections. Having determined the number of lineages using Maximum Likelihood and the species tree using coalescent analyses, we test their divergence time using a Bayesian approach. While the cytochrome b analysis only recovered two clades, the ddRAD analysis supported the reciprocal monophyly of five lineages of bald uakaris, with all clades including only individuals with distinct and exclusive diagnostic phenotypic characters. We found that species diversification in Cacajao occurred during the last 300 Kya and may have been influenced by the formation of rivers and flooded forests in western Amazonia. We propose that the four bald uakari subspecies currently recognised can be upgraded to species level and we describe the white uakaris from the basin of the Rio Tarauacá as a new species.

Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes), in Europe.

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.

Mapping differences in mammalian distributions and diversity using environmental DNA from rivers.

Finding more efficient ways to monitor and estimate the diversity of mammalian communities is a major step towards their management and conservation. Environmental DNA (eDNA) from river water has recently been shown to be a viable method for biomonitoring mammalian communities. Most of the studies to date have focused on the potential for eDNA to detect individual species, with little focus on describing patterns of community diversity and structure. Here, we first focus on the sampling effort required to reliably map the diversity and distribution of semi-aquatic and terrestrial mammals and allow inferences of community structure surrounding two rivers in southeastern England. Community diversity and composition was then assessed based on species richness and ß-diversity, with differences between communities partitioned into nestedness and turnover, and the sampling effort required to rapidly detect semi-aquatic and terrestrial species was evaluated based on species accumulation curves and occupancy modelling. eDNA metabarcoding detected 25 wild mammal species from five orders, representing the vast majority (82%) of the species expected in the area. The required sampling effort varied between orders, with common species (generally rodents, deer and lagomorphs) more readily detected, with carnivores detected less frequently. Measures of species richness differed between rivers (both overall and within each mammalian order) and patterns of ß-diversity revealed the importance of species replacement in sites within each river, against a pattern of species loss between the two rivers. eDNA metabarcoding demonstrated its capability to rapidly detect mammal species, allowing inferences of community composition that will better inform future sampling strategies for this Class. Importantly, this study highlights the potential use of eDNA data for investigating mammalian community dynamics over different spatial scales.

Fungal microbiomes are determined by host phylogeny and exhibit widespread associations with the bacterial microbiome.

Interactions between hosts and their resident microbial communities are a fundamental component of fitness for both agents. Though recent research has highlighted the importance of interactions between animals and their bacterial communities, comparative evidence for fungi is lacking, especially in natural populations. Using data from 49 species, we present novel evidence of strong covariation between fungal and bacterial communities across the host phylogeny, indicative of recruitment by hosts for specific suites of microbes. Using co-occurrence networks, we demonstrate marked variation across host taxonomy in patterns of covariation between bacterial and fungal abundances. Host phylogeny drives differences in the overall richness of bacterial and fungal communities, but the effect of diet on richness was only evident in the mammalian gut microbiome. Sample type, tissue storage and DNA extraction method also affected bacterial and fungal community composition, and future studies would benefit from standardized approaches to sample processing. Collectively these data indicate fungal microbiomes may play a key role in host fitness and suggest an urgent need to study multiple agents of the animal microbiome to accurately determine the strength and ecological significance of host-microbe interactions.

Space-time dynamics in monitoring neotropical fish communities using eDNA metabarcoding.

The biodiverse Neotropical ecoregion remains insufficiently assessed, poorly managed, and threatened by unregulated human activities. Novel, rapid and cost-effective DNA-based approaches are valuable to improve understanding of the biological communities and for biomonitoring in remote areas. Here, we evaluate the potential of environmental DNA (eDNA) metabarcoding for assessing the structure and distribution of fish communities by analysing water and sediment from 11 locations along the Jequitinhonha River catchment (Brazil). Each site was sampled twice, before and after a major rain event in a five-week period and fish diversity was estimated using high-throughput sequencing of 12S rRNA amplicons. In total, 252 Molecular Operational Taxonomic Units (MOTUs) and 34 fish species were recovered, including endemic, introduced, and previously unrecorded species for this basin. Spatio-temporal variation of eDNA from fish assemblages was observed and species richness was nearly twice as high before the major rain event compared to afterwards. Yet, peaks of diversity were primarily associated with only four of the locations. No correlation between ß-diversity and longitudinal distance or presence of dams was detected, but low species richness observed at sites located near dams might that these anthropogenic barriers may have an impact on local fish diversity. Unexpectedly high α-diversity levels recorded at the river mouth suggest that these sections should be further evaluated as putative "eDNA reservoirs" for rapid monitoring. By uncovering spatio-temporal changes, unrecorded biodiversity components, and putative anthropogenic impacts on fish assemblages, we further strengthen the potential of eDNA metabarcoding as a biomonitoring tool, especially in regions often neglected or difficult to access.

Fine-scale seascape genomics of an exploited marine species, the common cockle Cerastoderma edule, using a multimodelling approach.

Population dynamics of marine species that are sessile as adults are driven by oceanographic dispersal of larvae from spawning to nursery grounds. This is mediated by life-history traits such as the timing and frequency of spawning, larval behaviour and duration, and settlement success. Here, we use 1725 single nucleotide polymorphisms (SNPs) to study the fine-scale spatial genetic structure in the commercially important cockle species Cerastoderma edule and compare it to environmental variables and current-mediated larval dispersal within a modelling framework. Hydrodynamic modelling employing the NEMO Atlantic Margin Model (AMM15) was used to simulate larval transport and estimate connectivity between populations during spawning months (April-September), factoring in larval duration and interannual variability of ocean currents. Results at neutral loci reveal the existence of three separate genetic clusters (mean F ST = 0.021) within a relatively fine spatial scale in the north-west Atlantic. Environmental association analysis indicates that oceanographic currents and geographic proximity explain over 20% of the variance observed at neutral loci, while genetic variance (71%) at outlier loci was explained by sea surface temperature extremes. These results fill an important knowledge gap in the management of a commercially important and overexploited species, bringing us closer to understanding the role of larval dispersal in connecting populations at a fine geographic scale.

Environmental DNA metabarcoding as an effective and rapid tool for fish monitoring in canals.

We focus on a case study along an English canal comparing environmental DNA (eDNA) metabarcoding with two types of electrofishing techniques (wade-and-reach and boom-boat). In addition to corroborating data obtained by electrofishing, eDNA provided a wider snapshot of fish assemblages. Given the semi-lotic nature of canals, we encourage the use of eDNA as a fast and cost-effective tool to detect and monitor whole fish communities.

Environmental and anthropogenic drivers of connectivity patterns: A basis for prioritizing conservation efforts for threatened populations.

Ecosystem fragmentation and habitat loss have been the focus of landscape management due to restrictions on contemporary connectivity and dispersal of populations. Here, we used an individual approach to determine the drivers of genetic differentiation in caribou of the Canadian Rockies. We modelled the effects of isolation by distance, landscape resistance and predation risk and evaluated the consequences of individual migratory behaviour (seasonally migratory vs. sedentary) on gene flow in this threatened species. We applied distance-based and reciprocal causal modelling approaches, testing alternative hypotheses on the effects of geographic, topographic, environmental and local population-specific variables on genetic differentiation and relatedness among individuals. Overall, gene flow was restricted to neighbouring local populations, with spatial coordinates, local population size, groups and elevation explaining connectivity among individuals. Landscape resistance, geographic distances and predation risk were correlated with genetic distances, with correlations threefold higher for sedentary than for migratory caribou. As local caribou populations are increasingly isolated, our results indicate the need to address genetic connectivity, especially for populations with individuals displaying different migratory behaviours, whilst maintaining quality habitat both within and across the ranges of threatened populations.

Emerging Infectious Disease Implications of Invasive Mammalian Species: The Greater White-Toothed Shrew (Crocidura russula) Is Associated With a Novel Serovar of Pathogenic Leptospira in Ireland.

The greater white-toothed shrew (Crocidura russula) is an invasive mammalian species that was first recorded in Ireland in 2007. It currently occupies an area of approximately 7,600 km2 on the island. C. russula is normally distributed in Northern Africa and Western Europe, and was previously absent from the British Isles. Whilst invasive species can have dramatic and rapid impacts on faunal and floral communities, they may also be carriers of pathogens facilitating disease transmission in potentially naive populations. Pathogenic leptospires are endemic in Ireland and a significant cause of human and animal disease. From 18 trapped C. russula, 3 isolates of Leptospira were cultured. However, typing of these isolates by standard serological reference methods was negative, and suggested an, as yet, unidentified serovar. Sequence analysis of 16S ribosomal RNA and secY indicated that these novel isolates belong to Leptospira alstonii, a unique pathogenic species of which only 7 isolates have been described to date. Earlier isolations were limited geographically to China, Japan and Malaysia, and this leptospiral species had not previously been cultured from mammals. Restriction enzyme analysis (REA) further confirms the novelty of these strains since no similar patterns were observed with a reference database of leptospires. As with other pathogenic Leptospira species, these isolates contain lipL32 and do not grow in the presence of 8-azagunaine; however no evidence of disease was apparent after experimental infection of hamsters. These isolates are genetically related to L. alstonii but have a novel REA pattern; they represent a new serovar which we designate as serovar Room22. This study demonstrates that invasive mammalian species act as bridge vectors of novel zoonotic pathogens such as Leptospira.

Between the Balkans and the Baltic: Phylogeography of a Common Vole Mitochondrial DNA Lineage Limited to Central Europe.

The common vole (Microtus arvalis) has been a model species of small mammal for studying end-glacial colonization history. In the present study we expanded the sampling from central and eastern Europe, analyzing contemporary genetic structure to identify the role of a potential 'northern glacial refugium', i.e. a refugium at a higher latitude than the traditional Mediterranean refugia. Altogether we analyzed 786 cytochrome b (cytb) sequences (representing mitochondrial DNA; mtDNA) from the whole of Europe, adding 177 new sequences from central and eastern Europe, and we conducted analyses on eight microsatellite loci for 499 individuals (representing nuclear DNA) from central and eastern Europe, adding data on 311 new specimens. Our new data fill gaps in the vicinity of the Carpathian Mountains, the potential northern refugium, such that there is now dense sampling from the Balkans to the Baltic Sea. Here we present evidence that the Eastern mtDNA lineage of the common vole was present in the vicinity of this Carpathian refugium during the Last Glacial Maximum and the Younger Dryas. The Eastern lineage expanded from this refugium to the Baltic and shows low cytb nucleotide diversity in those most northerly parts of the distribution. Analyses of microsatellites revealed a similar pattern but also showed little differentiation between all of the populations sampled in central and eastern Europe.

Land-bridge calibration of molecular clocks and the post-glacial Colonization of Scandinavia by the Eurasian field vole Microtus agrestis.

Phylogeography interprets molecular genetic variation in a spatial and temporal context. Molecular clocks are frequently used to calibrate phylogeographic analyses, however there is mounting evidence that molecular rates decay over the relevant timescales. It is therefore essential that an appropriate rate is determined, consistent with the temporal scale of the specific analysis. This can be achieved by using temporally spaced data such as ancient DNA or by relating the divergence of lineages directly to contemporaneous external events of known time. Here we calibrate a Eurasian field vole (Microtus agrestis) mitochondrial genealogy from the well-established series of post-glacial geophysical changes that led to the formation of the Baltic Sea and the separation of the Scandinavian peninsula from the central European mainland. The field vole exhibits the common phylogeographic pattern of Scandinavian colonization from both the north and the south, however the southernmost of the two relevant lineages appears to have originated in situ on the Scandinavian peninsula, or possibly in the adjacent island of Zealand, around the close of the Younger Dryas. The mitochondrial substitution rate and the timescale for the genealogy are closely consistent with those obtained with a previous calibration, based on the separation of the British Isles from mainland Europe. However the result here is arguably more certain, given the level of confidence that can be placed in one of the central assumptions of the calibration, that field voles could not survive the last glaciation of the southern part of the Scandinavian peninsula. Furthermore, the similarity between the molecular clock rate estimated here and those obtained by sampling heterochronous (ancient) DNA (including that of a congeneric species) suggest that there is little disparity between the measured genetic divergence and the population divergence that is implicit in our land-bridge calibration.

The influence of habitat structure on genetic differentiation in red fox populations in north-eastern Poland.

The red fox (Vulpes vulpes) has the widest global distribution among terrestrial carnivore species, occupying most of the Northern Hemisphere in its native range. Because it carries diseases that can be transmitted to humans and domestic animals, it is important to gather information about their movements and dispersal in their natural habitat but it is difficult to do so at a broad scale with trapping and telemetry. In this study, we have described the genetic diversity and structure of red fox populations in six areas of north-eastern Poland, based on samples collected from 2002-2003. We tested 22 microsatellite loci isolated from the dog and the red fox genome to select a panel of nine polymorphic loci suitable for this study. Genetic differentiation between the six studied populations was low to moderate and analysis in Structure revealed a panmictic population in the region. Spatial autocorrelation among all individuals showed a pattern of decreasing relatedness with increasing distance and this was not significantly negative until 93 km, indicating a pattern of isolation-by-distance over a large area. However, there was no correlation between genetic distance and either Euclidean distance or least-cost path distance at the population level. There was a significant relationship between genetic distance and the proportion of large forests and water along the Euclidean distances. These types of habitats may influence dispersal paths taken by red foxes, which is useful information in terms of wildlife disease management.

Invading and expanding: range dynamics and ecological consequences of the greater white-toothed shrew (Crocidura russula) invasion in Ireland.

Establishing how invasive species impact upon pre-existing species is a fundamental question in ecology and conservation biology. The greater white-toothed shrew (Crocidura russula) is an invasive species in Ireland that was first recorded in 2007 and which, according to initial data, may be limiting the abundance/distribution of the pygmy shrew (Sorex minutus), previously Ireland's only shrew species. Because of these concerns, we undertook an intensive live-trapping survey (and used other data from live-trapping, sightings and bird of prey pellets/nest inspections collected between 2006 and 2013) to model the distribution and expansion of C. russula in Ireland and its impacts on Ireland's small mammal community. The main distribution range of C. russula was found to be approximately 7,600 km2 in 2013, with established outlier populations suggesting that the species is dispersing with human assistance within the island. The species is expanding rapidly for a small mammal, with a radial expansion rate of 5.5 km/yr overall (2008-2013), and independent estimates from live-trapping in 2012-2013 showing rates of 2.4-14.1 km/yr, 0.5-7.1 km/yr and 0-5.6 km/yr depending on the landscape features present. S. minutus is negatively associated with C. russula. S. minutus is completely absent at sites where C. russula is established and is only present at sites at the edge of and beyond the invasion range of C. russula. The speed of this invasion and the homogenous nature of the Irish landscape may mean that S. minutus has not had sufficient time to adapt to the sudden appearance of C. russula. This may mean the continued decline/disappearance of S. minutus as C. russula spreads throughout the island.

Preferred habitat and effective population size drive landscape genetic patterns in an endangered species.

Landscape genetics provides a framework for pinpointing environmental features that determine the important exchange of migrants among populations. These studies usually test the significance of environmental variables on gene flow, yet ignore one fundamental driver of genetic variation in small populations, effective population size, N(e). W(e) combined both approaches in evaluating genetic connectivity of a threatened ungulate, woodland caribou. We used least-cost paths to calculate matrices of resistance distance for landscape variables (preferred habitat, anthropogenic features and predation risk) and population-pairwise harmonic means of N(e), and correlated them with genetic distances, FST and D(c). Results showed that spatial configuration of preferred habitat and Ne were the two best predictors of genetic relationships. Additionally, controlling for the effect of Ne increased the strength of correlations of environmental variables with genetic distance, highlighting the significant underlying effect of Ne in modulating genetic drift and perceived spatial connectivity. We therefore have provided empirical support to emphasize preventing increased habitat loss and promoting population growth to ensure metapopulation viability.

A species-to-be? The genetic status and colonization history of the critically endangered Killarney shad.

Typically anadromous, the twaite shad (Alosa fallax) can become landlocked and adapt to a fully freshwater life. The only landlocked shad population in Northwestern Europe is found in a lake in Ireland, Lough Leane. The Killarney shad, Alosa killarnensis (or Alosa fallax killarnensis, as it is mostly referred to) displays a level of morphological divergence that indicates a long-term isolation in the lake. Microsatellites and mtDNA control region sequences were used within a coalescent framework (BEAST and Approximate Bayesian Computation (ABC)) to investigate its colonization history and to clarify its taxonomic status. Results indicate that the lake was likely colonized in two independent events, the first coinciding with the retreat of the ice sheet from the area after the Last Glacial Maximum and the second after the Younger Dryas. Microsatellite data has shown that these two landlocked lineages have completely admixed within the lake, and have started diverging significantly from their closest ancestor, the twaite shad. We argue that our molecular results (together with the life-history and physiological divergence between Killarney and twaite shad) are conspicuous enough to view the landlocked population as a new species, and one whose future existence would certainly not be insured by its sister taxon.

Human selection of elk behavioural traits in a landscape of fear.

Among agents of selection that shape phenotypic traits in animals, humans can cause more rapid changes than many natural factors. Studies have focused on human selection of morphological traits, but little is known about human selection of behavioural traits. By monitoring elk (Cervus elaphus) with satellite telemetry, we tested whether individuals harvested by hunters adopted less favourable behaviours than elk that survived the hunting season. Among 45 2-year-old males, harvested elk showed bolder behaviour, including higher movement rate and increased use of open areas, compared with surviving elk that showed less conspicuous behaviour. Personality clearly drove this pattern, given that inter-individual differences in movement rate were present before the onset of the hunting season. Elk that were harvested further increased their movement rate when the probability of encountering hunters was high (close to roads, flatter terrain, during the weekend), while elk that survived decreased movements and showed avoidance of open areas. Among 77 females (2-19 y.o.), personality traits were less evident and likely confounded by learning because females decreased their movement rate with increasing age. As with males, hunters typically harvested females with bold behavioural traits. Among less-experienced elk (2-9 y.o.), females that moved faster were harvested, while elk that moved slower and avoided open areas survived. Interestingly, movement rate decreased as age increased in those females that survived, but not in those that were eventually harvested. The latter clearly showed lower plasticity and adaptability to the local environment. All females older than 9 y.o. moved more slowly, avoided open areas and survived. Selection on behavioural traits is an important but often-ignored consequence of human exploitation of wild animals. Human hunting could evoke exploitation-induced evolutionary change, which, in turn, might oppose adaptive responses to natural and sexual selection.